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ABSTRACT 

We describe planetesimal accretion calculations for the dusty ring observed in the nearby AO star HR 
4796A. Models with initial masses of 10-20 times the minimum mass solar nebula produce a ring of 
width 7-15 AU and height 0.3-0.6 AU at 70 AU in 10 Myr. The ring has a radial optical depth ~ 1. 
These results agree with limits derived from infrared images and from the excess infrared luminosity. 

Subject headings: planetary systems - solar system: formation - stars: formation - stars: individual 
(HR 4796A) - circumstellar matter 



1. INTRODUCTION 

HR 4796A is a nearby A star with a large infrared (IR) 
excess. Jura (1991) measured the far-IR excess of this wide 
binary using IRAS data. Jura et al. (1995, 1998) associ- 
ated the excess with the AO primary and derived the ratio 
of the far-IR to stellar luminosity, Lpm/Li, « 5 x 10"^. 
In 1998, two groups reported extended thermal emission 
at A = 20 /im from a dusty disk with an inner hole at ~ 
40-70 AU assuming the Hipparcos distance of 67 ± 3.5 pc 



( Jayawardhana et al. 1998; Koerner et al. 1998). Obser 



vations with NICMOS aboard HST have revealed a thin 
annulus of scattered light, with a width of < 17 AU at a 
distance of ^ 70 AU from the centr al star ( Schneider ej 
al. 199SI). With an age of ~ 10 Myr ( ^tauffer al. 1995 



Barrado y Navascues et al. 1997), the AO star is older than 



most pre-main-sequence stars and younger than stars like 
/3 Pic and a Lyr with 'debris disks'. 

The dusty ring in HR 4796A challenges theories of 
planet formation. In most planetesimal accretion calcula- 
tions, planet-sized objects do not form on short timescales 
at large distances from the central star. Kenyon & Luu 
(1999; KL99 hereafter) estimate formation times of 10-40 
Myr for Pluto at 35 AU from the Sun. Achieving shorter 
timescales at 70 AU in HR 4796A requires large initial 
masses, which might conflict with masses derived from IR 
observations. In the inner Solar System, planet formation 
cannot be confined to a narrow ring, because high veloc- 
ity objects in adjacent ann uli interact and 'mix' planetar y 
growth over a large area ( Weidenschilling et al. 1997). 
This problem may be reduced at larger distances from the 
central star, where planetary growth is "calmer". 

Our goal in this paper is to develop planetesimal ac- 
cretion models that can lead to the dusty ring observed 
in HR 4796A. We begin in §2 with Monte Carlo calcu- 
lations to constrain the geometry and optical depth of 
dust in the ring. In §3, we derive plausible initial con- 
ditions which produce the observed dust distribution on 
10 Myr timescales. These models also satisfy constraints 



on the dust mass from IRAS observations and lead to a 
self-consistent picture for ring formation. We conclude in 
§4 with a brief summary and discussion of the implications 
of this study for planet formation in other star systems. 

2. MODEL IMAGES 

Current data constrain the geometry and optical depth 
of the ring. Near-IR images measure the amount of scat- 
tered light from the ratio of the 1.1-1.6 /im ra diation to 
the stellar luminosity, Lnir/L^, w 2 x 10~^ ( pchneidei 



et al. 1999). The far-IR luminosity limits the amount of 
stellar radiation absorbed and reradiatcd. To construct 
a physical model, we assume an annulus of width Aa 
and height z at a distance a = 70 AU from the central 
star. The luminosity ratios depend on the solid angle 
n/Air — 27raz/47ra^ — z/2a, the radial optical depth r, 
and the albedo u: Liqin/Li, = TU!{z/2a) and Lpiji/Li^ 
= t(1 — Lij){z/2a). These equations assume gray opacity 
and scattering in the geometric optics limit. If the annulus 
contains pla netesimals and dust in dynamical equilibrium, 
z/Aa <^ 1 (Hornung et al. 1985). Anticipating the re- 
sults of our coagulation calculations, where z/a ~ 10~^, 
we then have m k, Q.Z - clo se to observed values in (3 Pic 
(Packman & Paresce 1993) - and t ~ 1. 



We const ruct scattered light ima ges using a 3D Monte 
Carlo code (Wood & Reynolds 1999) with forced first scat- 
tering (Witt 1977) and a "peeling-off" procedure (Yusef- 
Zadeh, Morris, & White 1984). We adopt a dust num- 
ber density, n = nQe~^ g-(a-70) /2A ^ -^^j^gj-g ^j^g gcale 
height H and scale length A are in AU. We assume uj 



0.3, i sotropic scattering (see Figure 12 of Augereau et al. 



1999), and adjust r until the models yield Ljv//j/L* 
1.5 X lO^'^ for an input H and A. Model images with ojt 
— constant are identical in the optically thin limit. 

Figure 1 compares several models with the NICMOS 
1.1 /im image (from FITS data kindly sent by G. Schnei- 
der). We convolved Monte Carlo images with a gaussian 
point-spread function with FWHM = 0'.'12 to approximate 
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the 0'.'12 resolution of NICMOS ( Schneider et al. 19991 ). 
Model images with H > 5 AU (FWHM = 14 AU) or A > 
10 A U (FWHM ^ 27 AU) are more extended than the 
data ( Augereau et al. 1999| ). Our preferred model with 
H = 0.5 AU, A = 5 AU, and ujt = 0.25 reproduces the 
size and shape of the NICMOS image as well as the limb 
brightening observed towards the ring edges. These results 
match NICMOS flux ratios best for our adopted geometry; 
larger H implies smaller ujt. The 3(T Umit, ujt — 0.12-0.35 , 



agrees with previous estimates (cf. Koerner et al. 1998 



Schneider et al. 1999; Augereau et al. 1999| ). We disa gree, 
however, with the r ^ lO^'^ of Schneider et al.; their result 
is valid only for scattering in a spherical shell. 

3. COAGULATION MODEL 

To calculate dust evolution in HR 4976A, we use a 
coagul ation code based on the particle-in-a-box method 
(KL9£). This formalism treats planetesimals as a statis- 
tical ensemble of bodies with a distribution of horizontal 



and vertical velocities about Keplerian orbits (Safrono^ 
1969). We begin with a size distribution of Ni bodies hav- 
ing total mass Mi in each of i mass batches. Collisions 
among these bodies produce (i) growth through mergers 
along with cratering debris for low impact velocities or 
(ii) catastrophic disruption into numerous small fragments 
for high impact velocities. Inelastic collisions, long range 
gravitational interactions (dynamical friction and viscous 
stirring), and gas drag change the velocities of the mass 
batches with time. The code has been tested against an- 
alytic solutions of the coagulation equation and published 
calculations of planetesimal growth. Although inappropri- 
ate for the last stages of planet for mation, our approac h 
well-approximates the early stages (Kokubo & Ida 1996). 



We model planetesimal growth in an annulus of width 
Aa = 12 AU centered at a = 70 AU. The central star has 
a mass of 2.5 Mq. The input size distribution has equal 
mass in each of 38 mass batches with initial radii = 
1-80 m. For a Minimum Mass Solar Nebula with mass 
Mmmsn, the total mass in the annulus is Mq « 15 Me; 
the initial number of bodies with = 1 m is A^o ~ 3 x 10^'^. 
All batches start with the same initial velocity. The mass 
density po = 1.5 g cm~'^, intrinsic strength Sq = 2x 10^ erg 
g^^, and other bulk properties of the grains are adopted 
from earlier work (see KL99). 

Planetesimal growth at 70 AU follows the evolution 
described previously (KL99). The 80 m bodies first 
grow slowly into 1 km objects. During this slow growth 
phase, frequent collisions damp the velocity dispersion of 
all bodies. "Runaway growth" begins when the gravita- 
tional range of large objects exceeds their geometric cross- 
section. These bodies grow from 1 km up to ^ 100 km in 
several Myr. During runaway growth, coUisional debris, 
dynamical friction, and viscous stirring increase the veloc- 
ity dispersion of small bodies from 1 m s~^ up to 
40 m s~^. This evolution reduces the gravitational range 
of the 100 km objects and ends runaway growth. The 
largest objects then grow slowly to 1000-f km sizes. 

Figure 2(a) shows the growth of the largest object in 
several models. For Mq — 10 Mmmsn and cq = 10~^, 
Pluto-sized objects form in tp = 2.1 Myr at a = 35 AU, 13 
Myr at 70 AU, and 93 Myr at 140 AU. Models with smaller 
Mq take longer to make "Pluto" . Plutos form more quickly 
for Co < 10~^, because gravitational focusing factors are 



larger. 

Figure 2(b) shows the evolution of the scale height H 
for small objects. Initially, H — 27ra sin i < 0.003a for 
Co < 10~^. CoUisional damping cools the bodies during 
the slow-growth phase; H remains small. H increases dra- 
matically during runaway growth, when dynamical pro- 
cesses heat up the smallest bodies. Once runaway growth 
ends, H slowly increases to 0.3-0.6 AU independent of Mq, 
Co, and other input parameters. 

When H begins to increase, high velocity collisions pro- 
duce numerous "dust grains" with sizes ^ 1 m. We do 
not follow explicitly the evolution of these bodies. In- 
stead, we assume that coUisional debris is (i) swept up 
by 1 m or larger objects, (ii) ejected by radiation pressure, 
or (iii) dragged inwards by the Poynting- Robertson effect. 
Grains with sizes exceeding 4-5 /im are stable against ra - 
diation pressure ( Jura et al. 1998| ; Augereau et al. 1999). 
Poynting- Robertson drag reduces the mass in small grains 
on a timescale tpp « 1.0 Myr (ri/4 /im). With the short 
collision times, < 10^ yr, in our model annulus, 1 Myr 
seems a reasonable estimate of the timescale for collisions 
to produce 4 /im grains which are removed by radiative 
processes. For this paper, we calculate the accretion ex- 
plicitly and adopt a 1 Myr timescale for dust removal. 

Figure 2(c) shows the dust mass as a function of time. 
The results are not sensitive to the adopted mass distribu- 
tion for grains with ^ 4 ^m or to factor of 2-3 variations 
in the removal timescale. The dust mass is initially large 
due to the starting conditions. The dust mass decreases 
with time, because (i) coUisional damping of the smaller 
bodies leads to less coUisional debris and (ii) radiative pro- 
cesses and accretion by large bodies remove dust. Once 
runaway growth begins, collisions between small bodies 
produce more dust. The dust mass then reaches a rough 
equilibrium between collision debris and dust removed by 
radiation forces and by the larger bodies. 

These results indicate that large dust masses correlate 
with runaway growth and the formation of 1 or more Plu- 
tos in the outer parts of the disk. To predict the amount of 
radiation absorbed and scattered by dust and larger bod- 
ies, we compute t from the model size distribution. We 
assume the geometric optics limit because r; 3> A. For the 
large bodies t = X^i^i '^ifi^Q^i where rii is the number 
density in mass batch i, Ui is the extinction cross-section, 
and N is the number of mass batches. We adopt di — 27Trf 
and a volume Vi — 2TTaAaHi to compute rii — Ni/Vi and 
hence r for material with r.i ^ 1 m. 

Estimating t for small particles requires an adopted cu- 
mulative size distribution, Nc oc r~^. We consider three 
choices: (i) q = 2.5, the coUisional limit for coagulation; 
(ii) q — 3, equal mass per mass interval; and (iii) q = 3.5, 
the approximate distribution for grains in the interstellar 
medium. Our calculations produce q ~ 2.7 for 1-100 m 
bodies. We expect a slightly steeper mass distribution for 
smaller bodies, because collisions between smaller bodies 
produce fewer mergers and more debris. 

Figure 2(d) shows how r evolves for a single model. The 
large bodies initially have modest radial optical depth, 
tl ~ 0.2. This optical depth decreases with time, ex- 
cept for a brief period when runaway growth produces 
10-100 km objects with small scale height above the disk 
midplane. The large bodies are transparent once a Pluto 
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forms. The small grains are also initially opaque. This 
dust is transparent at late times if most of the mass is in 
the largest grains, q <^ 2.8. The dust is opaque for g <; 3. 

Table 1 summarizes results for various initial conditions. 
Models with Mq « 10-20 Mmmsn and eo « lO-^'-lO^^ 
achieve r ^ 1 in 10 Myr. Less massive disks produce less 
dust on longer timescales. The results are not sensitive 
to other input parameters, including the size distribution 
and the bulk properties of the bodies. 

Table 1 also shows why dust in HR 4796A lies in a ring. 
In disks with surface density E oc a~'^/^, the Pluto for- 
mation timescale^ is tp « 13 Myr (Mq/IO Mmmsn)^^ 
(a/70 AU)^'^. Once an annulus at a begins to form dust, 
material at a + Aa must wait a time, At/tp « 2.7Aa/a, to 
reach the same state. This result sets a hard outer limit to 
the ring, Aa/a w OAAt/tp « 0.1-0.2, if At is the time for 
H to double in size during runaway growth, 2-3 Myr. We 



expec ; a hard inner edge, because particle velocities reach 
the shattering limit of ^ 100 m s~^ (KL9 9) or planets 
sweep up the dust (e.g., Pollack et al. 1996) or both. 



4. DISCUSSION AND SUMMARY 



0\i \ results indicate that the dusty ring in HR 4796A is 
a natural outcome of planetesimal evolution. Planet for- 
mation at 70 AU in 10 Myr is possible with an initial disk 
mass of 10-20 Mmmsn- Dust production associated with 
planet formation is then confined to a ring with Aa ~ 
7-15 AU. The optical depth in this ring satisfies current 
constraints on scattered light at 1-2 /xm and on thermal 
emission at 10-100 /im if the size distribution of the dust 
is Nc oc r"'' with g ^ 3 for ^ 1 m. Models with disk 
masses smaller than 10 Mmmsn fail to produce planets 
and an observable dusty ring in 10 Myr. 

An uncertainty in our model is the timescale to pro- 
duce 1-80 m bodies from small dust grains in a turbulent, 
gaseous disk. Cuzzi et al. (1993) show that grains grow 
very rapidly once they decouple from eddies in the disk. 
The decoupling timescale depends on the unknown disk 
viscosity at 70 AU. 

Our model makes several observational predictions. We 
expect L]sfifi/Li, = constant for A < 5 /zm; current data 
are consistent with this prediction at the 1.5a level. Bet- 
ter measurements of the ring flux at A > 1.6 /xm would 



test our optical depth assumptions and yield interesting 
constraints on grain properties. Deep images at A > 10- 
20 /im with high spatial resolution should detect material 
outside the ring. We predict r « 0.1 in large bodies for 
a ^ 80 AU; the surface brightness and temperature of this 
material should decrease markedly with radius. This ma- 
terial should have negligible mass in small objects, because 
coagulation concentrates mass in the largest objects when 
H is small. We also expect a flux of dust grains into the 
central star, although we cannot yet compare quantita- 
tive predictions with observations. Future calculations of 
radiative processes within the ring will address this issue. 

Applying this HR 4796A model to other stars with cir- 
cumstellar disks is challenging due to small number statis- 
tics and unfavorable circumstances. Nearby companion 
stars probably influence th e dynamics of dus t y rings in 
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HD 98800 and HD 141569 ( Pirzkal et 

1999| ; [Lagrange et al. 2000). In HR 4796A, the M-type 



companion lies well outside the ring radius and cannot 
modify ring dynamics significantly. Older systems like /3 
Pic and a Lyr require time-dependent treat ment of dust 



to al low the ring to spread with time (e.g., Artymowicz 



1997). We plan to incorporate this time-dependent be- 



havior in future calculations to see whether the ring in 
HR 4796A can evolve into a debris disk (as in e.g. f3 Pic 
and a Lyr) on a timescale of 100-200 Myr. 

The main alternative to in situ ring formation at 70 
AU is migration of a planet formed at a smaller radius. 
Weidenschilling & Marzari (1996) show that gravitational 
interactions can scatter large objects into the outer disk 
in less than 1 Myr. Migration reduces the required ring 
mass by a factor of 10-100. However, the scattered body 
has a large eccentricity, e ~ 0.5. Dynamical friction might 
circularize the orbit in 10 Myr, but would induce large ec- 
centricities in smaller bodies. The width of the dusty ring 
would probably exceed observational constraints. Future 
calculations can address these issues. 

We thank B. Bromley for helping us run our code on the 
HP Exemplar "Neptune" at JPL and for a generous allot- 
ment of computer time through funding from the NASA 
Offices of Mission to Planet Earth, Aeronautics, and Space 
Science. 
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^Pluto is a handy reference: 1000+ km objects form roughly in the middle of the rapid increase in H which produces large dust masses. 
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Fig. 1. — Comparison of model scattered light images with HST data of HR 4796A. (a) upper loft panel: NICMOS coronagraphic image at 
1.1 ^m. (b) upper right panel: Model scattered light image with z = 0.5 AU, R = 5 AU, and wri'^ifi = 0.25. (c) lower left panel: As in (b) 
for ^ = 5 AU, ii = 10 AU, and ujtmir = 0.02. (d) lower right panel; As in (c) lor z = 1 AU, -R = 20 AU, and LOTf^m = 0.1. 
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Fig. 2. — Results for HR 4796A coagulation models, (a) lower left panel: maximum radius for models with eo = 10~^ at a = 35, 70, and 
140 AU; initial masses are listed in the legend, (b) upper left panel: scale height of small bodies for 70 AU models with eo and Mo as listed, 
(c) upper right panel: dust mass for 70 AU models with eo and Mo as listed, (d) lower right panel: optical depth for a 70 AU model with 
eo = 10-3 and Mq = IOMmmsn- 



Table 1 

Results of Planetesimal Accretion Calculations 



a 


Mo 


eo 


tp 


log Ts 


log Tl 


35 


100 


10-3 


2.1 


0.09 


-3.43 


70 


15 


10-* 


81.1 


-1.90 


-3.92 


70 


45 


10-* 


20.4 


-0.96 


-3.27 


70 


150 


10-* 


5.7 


-0.04 


-2.69 


70 


15 


10-3 


156.4 


-2.21 


-4.24 


70 


45 


10-3 


50.0 


-1.42 


-3.65 


70 


75 


10-3 


29.9 


-0.94 


-3.48 


70 


150 


10-3 


13.0 


-0.34 


-3.25 


70 


300 


10-3 


6.6 


-0.05 


-2.62 


140 


200 


10-3 


92.6 


-1.20 


-2.95 



Note. — a is the distance of the annulus from the central star in AU; Mq is the initial mass in 
the annulus in Me\ is the initial eccentricity of each mass batch; tp is the timescale in Myr to 
produce Pluto-sizcd objects; ts is the optical depth in dust when the first Pluto forms, assuming 
equal mass in dust per decade in radius; tl is the optical depth of the large bodies when the first 
Pluto forms. 



